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Abstract 


(“a  lattice  theory  is  developed  for  composite  materials  in  which 
the  matrix  and  inclusions  of  the  composite  are  idealized  as  a  gridwork 
of  thermoelastic  bars  attached  to  small  rigid  masses.  This  model 
generalizes \the  spring  and  mass  model  prominent  in  the  classical 
lattice  theiry  of  solids.  Replacement  of  the  springs by  elastic 
bars  provides  extra  internal  degrees  of  freedom  beyond  the  central 
force  model  of  the  lattice  theory.  The  lattice  representation  leads 
to  a  system  of  finite  difference  equations  of  motion  which  can  be  used 
to  solve  dynamic  and  thermal  stress  problems  numerically.  Equations 
for  stress  and  couple  stress  are  obtained,  which  are  identical  in  form 


to  those  of  the  linear 


theory  of  micropolar  thermoelasticity.  \  . 


iii 


Introduction 


Consider  a  matrix  containing  inclusions  which  are  of  nearly  equal 
size  and  which  are  spaced  relatively  uniformly  (spacing  L) ,  as  shown 
in  Figure  la.  Such  a  composite  material  can  be  represented  by  the 
lattice  model  shown  in  Figure  lb,  in  which  the  matrix  and  inclusions 
have  been  replaced  by  a  gridwork  of  massless,  flexible  bars  attached 
to  small,  rigid  bodies.  The  model  so  introduced  represents  a 
generalization  of  the  mass-spring  model  well  known  in  the  theory  of 
crystal  lattice.  By  introducing  elastic  bars  to  represent  the 
interparticle  forces  one  provides  extra  internal  degrees  of  freedom  for 
non-central  forces  and  couples.  The  model  is  therefore  expected  to 
encompass  certain  physical  phenomena  not  present  in  the  classical  theory. 
Indeed  this  turns  out  to.  be  the  case. 

In  order  that  the  lattice  model  represents  the  actual  composite 

material  in  some  optimum  sense,  "effective”  values  must  be  selected  for 

the  masses  and  the  mass  moments  of  inertia  of  the  bodies,  the  stiffness 

and  the  thermal  expansion  coefficients  of  the  bars,  etc.  For  example, 

the  effective  mass  of  each  rigid  body  might  be  taken  as  the  sum  of  the 

3 

mass  of  a  single  inclusion  and  the  mass  of  a  volume  L  of  the  matrix. 
Similarly,  effective  stiffness  coefficients  for  each  bar  could  be  found 
by  considering  the  volume  of  material  between  each  inclusion  and  the 
elastic  moduli  of  the  composite  (determined,  for  example,  using  the  classic 
rule  of  mixtures).  It  is  not,  however,  the  purpose  of  this  paper  to 
discuss  the  procedure  for  selecting  the  parameters  of  the  lattice  model. 
Rather,  the  aim  here  is  to  develop  the  equations  which  govern  the  deformation 
of  such  materials. 
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EQUATIONS  OF  MOTION 


A  repeating  unit  of  a  two-dimensional  lattice  model  is  shown  in 
Figure  2.  Interactions  between  nearest  and  next  nearest  neighbors  are 
transmitted  by  means  of  thermoelastic  bars.  The  bars  are  capable  of 
supporting  axial,  shear  and  bending  loads,  as  well  as  transmitting  heat. 
The  equations  of  motion  for  the  rigid  body  at  point  (jl,m)  are 


..  jl,m  _  Jl,m  f  T 
m  u  ^  =  F  (a  =  1,2) 

a  a 

r  i.m  _  ^  2.  ym 
j  4.3  =  M3 


(1) 


where  u  and  are  the  translations  and  rotation  of  the  lattice 

a  ^  m  i  m 

pbint  (Jl,m),  and  F^  ’  and  *  refer  to  the  resultant  forces  and  moment 

exerted  on  the  particle  by  the  adjoining  beam  elements.  The  constants 

m  and  J  denote  the  mass  and  mass  moment  of  inertia  of  the  body.  These 

forces  and  moment  result  from  both  temperature  changes  0  in  the  beams  and 

from  relative  displacements  (translations  and  rotations)  of  the  end  points 

of  the  beams.  The  forces  which  each  element  exerts  on  mass  m  may  be 

calculated*  from  the  stiffness-matrix  and  the  thermal  expansion  properties 

of  that  particular  element.  It  is  assumed  that  all  bars  connecting  nearest 

lattice  points  have  identical  properties,  and  likewise  all  bars  connecting 

next-nearest  points  are  identical.  Summing  the  contributions  of  all  eight 

it  m  m 

bars  gives  the  resultant  generalized  forces  F^  ^  ^  .  Substitution 

of  these  expressions  into  (1)  yields 


m  u. 


iyva 


^  ,  £,+l,m  „  ,  i-lym. 

La^  (u^  -2u^  ’  +  u^  ^  ) 

,  ,  jl,m+l  ^  ,  Zym-1. 

+  Lh^  (u^  ’  -2u^  '^l  ^ 


*The  details  of  computing  the  forces  which  a  typical  bar  exerts  on  the  mass 
are  given  in  the  Appendix. 
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+  -2  (^2  +  ^2^  ^^1  ^1  "^1  -  4u^  ’  ) 


,  L  ^  -L  \  /  il+l,ni+l  .  jl+l,iii-l  £-l,m+l 
+  -  (a^-  b2)  (U2  -  ^^2  -  ^^2  +  ^^2  > 


+  —  bi  (<1.3 


£,in+l 


,  ,  ,,  Jl+l ,111+1  £+l,m-l  .  ,  iZ/-l,m+l  £-l,m-lv 

*3  “  *3  ^ 


Li 

2 


2 


g2 


^^^,+1,111+1  _j_  Qj?,+l,in-l  _  Qil“l,m+1  __  Qil-l,in-l^ 


m  u. 


il,ni 


_  .  jl,in+l  «  ,  £,in-lv 

=  La^  (u^  -  ’  ) 

.  Ti.  /  £.+l,m  «  jl,in  ,  Jl-l,niv 

+  Lb^  (u^  "  ’  ) 


,  L  .  .  .  X  /  Jl+l,Ta+l  ,  il+l,m-l  ,  l-lym+1  ^  il-l,m-l  , 

+  Y  (a2  +  b^)  (u^  +  -  4u2  M 


,  L  ,  u  il+l,m+l  j?,+l,ni-l  il-l,m+l  ,  il“l,tn-lv 

+  ^  (a2  -  b^)  "  ^1  “  ’  ) 


•.  /.  Jt+l,m  .  £-l,m. 

"  F  ^<^3  ■  ‘•’3  ^ 


L  ,  ,,  jl+l,m+l  £+l,m-"l  £— l,tn+l  ,  £— l,in—lv 

~  02(4)3  ^  ^3  ”  S  ”  ^3  ^ 
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_  Jki 

2  H  ^ 


-)l+l,in-l  „S,-l,m+l 

—  g2(0  -0  +0  -0  ) 


...(2) 


J<f3  =  — [(U2  ’  -  U2  ’  )  -  (u^  ’  -  ’  )] 


+  ,  r  .  JH-l,in+l  . 

-b2[(u2  +«2  -’^2  -'^2  ) 


(u^  -  ’^1  -  ’  )] 


-  I  +  4.3^"^’"*  +  +  43''’”‘^  +  8<|>3"*”] 


"  7  ^2  *-^3  ^*^3  +4^3  +  8(1)^  ] 


Here  the  stiffness  coefficients  a  ,  b  ,  d  ,  are  defined  as 
E  A 

a  - 

5  ■  L-Lc 


12E  I 

b  =  - LJ. 

?  -  L.L^3 


5  =  1,2 


(no  summation) 


(3) 


d.  E 


?  -  L*L 


n  ■  L-L^ 
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where  E^,  A^,  a  ,  and  L  denote  the  Young's  modulus,  cross-sectional 
area,  moment  of  inertia,  coefficient  of  thermal  expansion,  and  the  length 
of  a  bar  respectively.  A  subscript  ^  =  1  designates  bars  connecting 
nearest  lattice  points;  5=2  refers  to  elements  connecting  next -nearest 
points. 

The  equations  of  motion  (2)  can  be  written  as  differential  equations 
by  first  expanding  the  displacements  of  neighboring  lattice  points  in 
Taylor  series  about  the  point  (jt,m);  i.e., 

u  =  (1  +  L3.  +  ^  9  2  +  ...)  u 

a  1  2!  1  a 


u 


il-l,m 


l2 


=  (1  -  L9^  +  2r  -  •••>  "a 


£,m 


(4) 


u  =  (1  +  l9,  +  L9„  +  ^  9,2  +  l29-9„  +  ^  9,2  +  ...)  u  *■’ 

a  12  2!  1  12  2!  2  a 


m 


u  =  (1  +  L9^  -  L9„  +  ^  92  _  L2a  3  +  ^  9,2  +  ...)  u 

a  1  2211  12212  a 


m 


u  =  (1  -  L9i  -  L9,  +  —  9/  +  L^9,9,  + 


^  ~2j.T22'i  j.  ii —  g  2  _ 


1“2  ■  2!  “2 


)  u 


il,m 


where  9^  and  9^  respectively  represent  the  partial  differential  operators 
3/9x^  and  9/9x2- 

Substituting  (4)  along  with  similar  expansions  for  the  rotations 
into  (2),  and  retaining  only  second  derivatives,  gives 


^  Ui  =  (a^  +  a2  +  h^)  +  b2)  u^^22 


+  (2a2  -  2b ^2,12'^  +  ^^2^  ‘^3,2 


(s^  "b  2^2^09^ 


(5) 


5 


m 


u. 


(a^  +  a 
+  (2a2 


2  ^2^  ^2,22 


■“  ^^2^  ^1,12 


+  (32  +  +  b2)  U2^ 

•  +  2b2)<l>3^ 


-  (gj^  +  2g2)0»2 


J_ 


(bi  +  2b2)(u2^^  -  u^^2> 


+  d2)  ((|) 


3,11 


'3,22^  " 


(5) 


Equations  (2)  and  (5)  provide  a  set  of  difference  and  differential 
equations,  respectively,  which  govern  the  motion  of  the  two-dimensional 
lattice  model  under  consideration.  A  solution  to  these  equations  which 
satisfies  prescribed  boundary  conditions  yields  the  displacements  u^ 
and  the  rotation  cf)^.  One  can  then  compute  the  forces  and  moments  in  the 
lattice  bars  using  the  appropriate  stiffness  coefficients  (given  in  the 
Appendix).  Finally,  stresses  and  couple  stresses  can  be  found  by 
averaging  the  forces  and  moments  over  the  faces  of  a  small  cube  surrounding 
a  lattice  point. 


STRESSES  AND  COUPLE  STRESSES 


To  obtain  the  general  constitutive  equation  for  the  material 
represented  by  this  model  we  isolate  the  cube  surrounding  mass  (^,m), 
as  shown  in  Fig.  3a.  This  element  is  subject  to  the  concentrated  forces 
and  moments  transmitted  by  those  bars  which  are  intersected  by  the  cube. 
These  generalized  forces  are  expressed  in  terms  of  the  translations  and 
rotations  at  the  end-points  of  the  beams,  as  described  in  the  Appendix. 
Averaging  the  forces  over  the  faces  of  the  cube,  and  then  making  use  of 
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(4)  where  only  first  derivatives  are  retained,  gives  the  stress  components 
(see  Fig.  3 


^11  °  ^2  ^  ^2^  '^1,1  ^^2  "  ^2^  '^2,2  ■  ^®2^® 

^22  °  ^2^  '"2,2  ^^2  "  ^"2^  ""l,!  "  ^  ^®2^® 

’^12  =  S  ■  ^2^  ""1,2  ^^2  +  ’’i  +  ’’2^  “2,1  “  ^’"l 

^21  °  ^^2  "  ^2^  “2,1  ^^2  ’"l  ^2^  “l,2 

and  the  couple-stress  components 
di  d2 

”*13  =  +  ~^*3,1 

”23  =  +  T^‘'’3,2 


(6) 


(7) 


It  is  of  interest  to  compare  equations  (6)  and  (7)  with  the 
constitutive  equations  of  Eringen’s  micropolar  theory  [1,2].  The  stresses 
and  couple  stresses  in  a  linear  micropolar  anisotro^j^  thermoelastic 
solid  are  given  by 


^jlkmn  ^mn 
™k£  ””  ^£kmn  ^m,n 


(8) 


where 


^jlkmn  \in£,k’  ^jtkmn  ^mn^k 


(9) 
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and 


\)l  -  2 

1 

r  E  OB  u 
m  ^  mnp  p,n 


(10) 


In  equations  (8)  and  (10)  (fi  is  the  "microrotation", 

m 

rotation  of  "microelements”  about  the  center  of  mass 
"macrovolume  element".  For  the  present  model  we  can 
rigid  body  at  lattice  point  (£,m)  as  a  microelement, 
cube  as  a  macrovolume  element. 

In  view  of  (9)  the  number  of  coefficients 

itkmn 

For  an  orthotropic  material  (one  having  3  orthogonal 
equations  (8)  become,  in  the  two-dimensional  case 


i.e.,  the  average 
of  a  deformed 
interpret  the  small 
and  the  surrounding 

planes  of  symmetry) 


^11  "^1111  ^'ll  ^1122^22  ®11® 

*^22  "  ^1122^11  ^2222^22  ®22® 

hi  "  ^2112^12.  ^2121^21 

^21  °  ^1212^12  ^2112^21 


(11) 


and 


”*13  "  Sl3l't’3,l 


(12) 


”*23  S232‘<’3,2 


The  lattice  model  under  consideration  is  symmetric  also  with  respect  to 
the  two  orthogonal  45®  planes  which  bisect  the  planes  of  symmetry 
and  x^  -  x^ 


8 


(see  Fig.  2).  Accounting  for  this  symmetry  and  using  (10),  equations 
(11)  reduce  to 


^11  °  ^1122'^2,2  ®11® 

^22  °  ■^1122'^!, 1  ^llll'^2,2 

^12  "  ■^2112“l,2  ^1212“2,1  "  ^^1212  “  ^2112^ ‘•’3 

^21  "  ■^1212“l,2  ^2112'^2,1  ^^1212  "  ^2112^ ‘•’s 

and  equations  (12)  become 


“l3  "  ^313l‘*’3,l 
”*23  “  Sl3l'*’3,l 


(14) 


A  comparison  of  (6),  (7)  with  (13),  (14)  shows  that  the  constitutive  equations 
for  the  lattice  model  and  the  micropolar  theory  are  identical  providing 


^1111  =  +  ^2  ^2 

\l22  ^  ^2  ~  ^2 

^2112  "  ^2  "  ^2 

^1212  °  ^2  ’^l  ^2 

=  -  (g^  +  2g2) 

^1  *^2 

r  ~  — =.  j _ ^ 

3131  42 


(15) 
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Hence,  the  material  constants  in  the  micropolar  theory  can  be  evaluated 
in  terms  of  the  stiffness  and  thermal  expansion  properties  of  the  lattice 
elements.  Considerable  insight  into  the  relative  magnitudes  of  the 
material  constants  is  thus  provided.  The  approach  also  suggests  a 
method  for  investigating  how  the  elastic  properties  of  a  micropolar 
solid  would  be  affected  by  changes  in  the  geometrical  and  material 
properties  of  the  constituents.  In  effect,  it  permits  the  construction 
of  a  material  with  a  desired  property,  within  the  context  of  micropolar 
theory. 

The  equations  of  motion  for  the  element  of  mass  (i.e.,  the  cube) 
are  given  by  the  first  and  second  laws  of  motion  of  Cauchy  which  express 
the  local  balance  of  momenta. 

Restricted  to  the  two-dimensional  case,  in  the  absence  of  body  forces 
and  couples  these  are 


and 


pii 

a 


(ct,6  =  1,2) 


(16) 


“33,3  ^33y‘^3y  ^^‘*’3 


(17) 


Substituting  (6)  into  (16)  gives  equations  (5^  2^  providing 


m 

P  =  7T 


(18) 


Hence  the  displacement  equations  of  motion  for  the  lattice  point 
(^,m)  are  identical  to  those  for  the  cube  element. 

Substituting  (6)  and  (7)  into  (18)  gives 


Pj<l>3  -  (b^  + 

^1  ^2 

^”■'■“^^'^3,11  '*’3,22^ 


(19) 


-  (2b^  +  4b2)((i3 
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Letting 


Pj  ^  (20) 

we  note  that  equation  (19)  then  differs  from  (5^)  only  in  the 
coefficient  of  the  term  involving  (cf)^  difference 

is  due  to  the  fact  that  the  concentrated  couples  actii^g  on  the  lattice 
point  (i^,m)  are  different  from  those  acting  on  the  cube  element.  Note 
that  the  cube  is  intersected  by  some  lattice  bars  (e.g.  bars  (X,,m+1  - 
)l+l,m),  (Jl+l,m  -  il,m--l),  etc.)  which  do  not  intersect  the  lattice 
point.  These  bars  therefore  transmit  couples  to  the  cube,  but  not  to 
the  lattice  point.  In  addition  there  are  couples  acting  on  the  cube 
which  are  associated  with  the  transfer  of  shear  forces  from  the  lattice 
point  to  the  faces  of  the  cube. 
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APPENDIX;  INTERACTION  FORCES  IN  THE  LATTICE  MODEL 


The  nature  of  the  interaction  forces  transmitted  by  a  thermoelastic 
beam  in  the  lattice  model  is  now  examined.  Each  element  is  assumed  to 
be  a  straight  bar  of  uniform  cross-sectional  area  A,  flexural  rigidity 
El,  length  L  and  coefficient  of  thermal  expansion  a.  An  element  is 
capable  of  resisting  an  axial  force,  shear  and  flexure.  The  location  and 
direction  of  the  generalized  forces  (Q^)  the  corresponding  generalized 
displacements  (q^)  which  act  on  a  typical  element  are  shown  in  Fig.  4. 

The  forces  are  linearly  related  to  the  displacements  and  to  the 
average  temperature  change  A0  in  the  beam.  The  temperature  is  assumed  to 
be  constant  over  any  cross-section,  and  to  vary  linearly  along  the  axis 
of  the  beam.  In  this  case  the  forces  are  given  by* 

{Q}  =  [S]  {q}  +  [S^]A0  (A.l) 


*A  rectangular  matrix  is  denoted  by  [  ];  a  column  matrix  by  {  }. 
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and  the  forces  resulting  from  thermal  expansion  of  the  bars  are 
determined  using 


[S^] 


aEA 

0 

0 

aEA 


0 

0 


(A.  3) 


We  now  consider  a  particular  beam  element,  for  example  the  one  connecting 
lattice  points  (il,m)  and  +  1,  m  +  1)  as  shown  in  Figure  5,  P^, 

P^  and  P^  denote  the  generalized  forces  which  act  on  mass  (Jt,m),  and 
P^,  P^  and  P^  denote  those  acting  on  particle  (£  +  1,  m  +  1).  These 
forces  may  be  expressed  in  terms  of  the  equal  and  opposite  forces  which 
the  particle  exerts  on  the  beam 


{P}  =  [A]{Q} 


where  from  geometry 


ri  fi 

1  _  JL 

ri  ' 


0  0  -  1 

0  0  0 

0  0  0 

0  0  0 


(A.4) 


0 


0 


0 

1 

ri 

■  1 

0 


0 


0 


0 

1 

fl 

1 

fl 

0 


0 


0 


0 

0 


0 


-  1 


(A.4’) 


\ 


y 

I 
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The  generalized  displacements  (assumed  infinitesimal)  may  be 
expressed  in  terms  of  the  displacements  of  the  masses  (£ ,m)  and 
(Jl  +  1,  m  +  1)  as 


{q}  =  [B]{u, 


il,m  il,m 


H  ,m 

’3  ’ 


u. 


)l+l,m+l  £+l,m+l 


u, 


,  iH-l,m+L 

<^3  ^ 


(A.5) 


where 


[B] 


1 

1 

ri 

0 

0 

0 

0 


1 

n 

1 

ri 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

1 

1 

\J 

fl 

ri 

0 

_ 

fl 

1 

n 

0 

0 

0 

0 

0 

0 

0 

0 

1 


(A.5*) 


The  average  temperature  change  in  the  beam  is 


A0  - 


[cHe*'’”, 


Qjl+l,in+lj 


(A.  6) 


where 


[C]  =  [1/2  1/2] 


(A.6') 
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Substituting  (A.  2)  through  (A. 6)  into  (A.l)  gives 


r-nT  rxrir  ±  il+l,m+l  il+l,in+l  , 

iP}  =  [K]{u^  ’  ,  u^  ’  ,  <|>3  ’  ,  u^  ,  u^  ,  (J)^  } 


(A.  7) 


where 


[K]  =  [A][S][B] 
[K^]  =  [A][S^][C] 


(A.7’) 


In  expanded  motion,  the  interaction  forces  transmitted  to  mass 
(il,m)  become 

L  .  .  V  \/  il+l,m+l  jl,m. 

?!  =  -  (a2  +  b2)(u^  -  “l  > 


,  L  ,  V  \  /  ^+1  ,m+l  i  ,m. 

+  -2  (^2  -  (u^  "  ^^2  ^ 


+  f-  b2  (<^3 


il+l,m+l  _j_  ^  Jl,m^ 


^  g2 


(A.  8) 


_  L  r  I  -u  \ /TT  i^+l»ni4-l 

P2  =  -2  (^2  ^2^  ^^2  “  ^2  ^ 


,  L  /  n  N  /  il+l,m+l  Jl,m. 
+  -(a2-b2)(u^  -Ui’) 


-  ^  b2  +  <^3^’”) 


f  g2  (0^+1 ’“+1  +  0^’"*) 
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P3  -  —  b2  [(U2  -  U2  )  -  (u^  -  ’  )] 


i  ^  /•*  i!-+l,m+l  „  £,m. 


where 


= 


E2A2 


2  L-L, 


’’2  = 


L-L23 

“2^2 


‘2  L-L. 


(A.  9) 


L-L. 


The  subscripts  2  in  (A.  8)  and  (A.  9)  have  been  introduced  in  order 
to  indicate  that  these  quantities  refer  to  a  beam  element  which  connects 
next-nearest  lattice  points.  A  subscript  1  is  used  to  designate  nearest 
lattice  points. 

Expressions  for  the  interaction  forces  in  the  other  bars  are  found 
in  a  similar  manner. 
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